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DETERMINATION OF TIME DEPENDENT FACTORS OF 
COEFFICIENTS IN FRACTIONAL DIFFUSION EQUATIONS 


KENICHI FUJISHIRO * 1 ) AND YAVAR KIAN 2 ) 

Abstract. In the present paper, we consider initial-boundary value problems for partial 
differential equations with time-fractional derivatives which evolve in Q = SI x (0, T) where 
Pis a bounded domain of R d and T > 0. We study the stability of the inverse problem of 
determining the time-dependent parameter in a source term or a coefficient of zero-th order 
term from observations of the solution at a point Xq £ SI for all t £ (0 ,T). 


1. Introduction 

Let Q be a bounded domain of M. d , d — 1,2, 3, with C 2 boundary <9fi. We set £ = dQx (0, T ) 
and Q = x (0, T). We consider the following two initial-boundary value problem (IBVP 
in short) for the fractional diffusion equation 

{ d^u(x, t ) + Au(x, t) = f(t)R(x , t), (x, t) G Q, 

B a u(x,t) — 0, (x,t)e£, (1.1) 

u(x, 0) = 0, x G 

and 

{ d^v(x, t ) + An (a:, i) + f(t)q(x, t)v(x, t ) = 0, (x, t) G Q, 

B a v(x,t ) = 0, (x,t) G £, (1.2) 

v(x, 0) = v 0 (x), x G 

with 0 < a < 1. Here we denote by <9“ the Caputo fractional derivative with respect to t: 

5 ° s(() ;= r(T^)/ ( ‘- s) "”5f (s)ds - 

The differentiual operator A is defined by 

Au(x,t) := - > 

i,j=l 1 ' ' 

where the coefficients satisfy 

d 

a ij = a,ji, 1 < i,j < d, and a ij(x)£i£j > /t|£| 2 , x G ff, £ G IR rf 
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for some (i > 0. Moreover B a is defined as 

B a u(x) = (1 — a(x))u(x ) + a(x)d UA u(x), x G d fi, 


where 


d UA u(x) = 

i,j =1 


ai A x )^A x ) 


and v = (yii ..*% I'd) is the outward normal unit vector to <912. Here o is a C 2 function on <912 
satisfying 

0 < cr(x) <1, x G <912. 

For the regularity of a^-, we assume 


dij G £9(12) if a = 0, 
a,ij G C 2 (12) if cr ^ 0. 


Note that the regularity for a tJ depends on whether a = 0 or not, which is due to condition 
(2.3) in the next section. 

In the present paper, we consider the inverse problem which consists of determining the 
function {f (t)} t e(o,T) hi (1-1) and (1.2) from the observation of the solution at a point x 0 G 12 
for all t G (0, T). 

The partial differential equations with time fractional derivatives such as (1.1) and (1.2) are 
proposed as new models describing the anomalous diffusion phenomena. Adams and Gclhar 
[1] pointed out that the held data in a highly heterogeneous aquifer cannot be described 
well by the classical advection diffusion equation. Hatano and Hatano [11] applied the 
continuous-time random walk (CTRW) as a microscopic model of the diffusion of ions in 
heterogeneous media. From the CTRW model, one can derive a fractional diffusion equation 
as a macroscopic model (see e.g., Metzler and Klafter [15] and Roman and Alemany [18]). 
In particular, the fractional diffusion equation can be used as a model for the diffusion of 
contaminants in a soil. Therefore the inverse problem considered in this paper means the 
determination of the time evolution of pollution source in (1.1) and reaction rate of pollutants 
in (1.2) respectively. In this paper, we consider such problems assuming the boundedness of 
the time-dependent parameter {/(t)} tG ( 0 ,r) (see (2.1)). 

As monographs of fractional calculus, there are books such as Podlubny [16] and Samko, 
Kilbas and Marichev [22] for example. As for mathematical works concerned with partial 
differential equations with time fractional derivatives, we can refer to Agarwal [3], Gejji and 
Jafari [9], Gorenflo and Mainardi [10], Luchko [14] and references therein. 


The remainder of this paper is composed of four sections. In Section 2, we state our main 
results. In Section 3, we study the forward problem for the IB VPs (1.1) and (1.2) and prove 
the unique existence and regularity of the solutions. In Sections 4 and 5, we complete the 
proof of our main results—Theorems 2.1 and 2.2 respectively. 
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2. Main results 

By L 2 (Q), we denote the usual L 2 -spaee equipped with the inner product (•, •) and the 
norm || • || := || • Moreover H S (Q), seR, and W m ' p (Ll), m G N, 1 < p < oo, are the 

Sobolev spaces (see Adams [2] for example). 

For the time dependent parameter {f(t)}te( o,t), we always assume 

/ G L°°(0,T). (2.1) 

For other given functions in (1.1), we suppose 

R G L p ( 0, T; H 2 (n)), - < p < oo and B a R = 0 on E. (2.2) 

a 

On the other hand, in the IB VP (1.2), we suppose 

fg G L°°(0, T; H 2 (Q)) (and d v q = 0 on E if a ^ 0), 

|uo G and B a v o = B a (Av o) = 0 on <90. 

Assuming these conditions, we prove in Section 3 that the IB VPs (1.1) and (1.2) admit unique 
solutions it, v G C([0, T]; i/ 2 (0)) with d/u G L p (0, T; H S (Q)) and G L°°(0, T; H S (Q)) for 
some s > d/2. Therefore, using the Sobolev embedding theorem (see Theorem 9.8 in Chapter 
1 of [13] for example), for any xq G O, we see that 

d/u(x o, •) G L p (0, T) and d?v(x 0 , ■) G L°°(0, T). 

Then our main results can be stated as follows; 


Theorem 2.1. Let condition (2.2) be fulfilled and Ui be the solution of (1.1) for f — f t G 
L°°(0,T), (i = 1,2). We assume that there exist x$ G 0 and 5 > 0 such that 

\R(x 0 ,t)\>5, a.e. 1 G (0,T). (2.4) 

Then there exists a constant C > 0 depending onp, T, 0, 5 and ||-R||LP(o,T ; ir 2 ( 0 )) such that 

ll/i — /2||lp(o,t) < C\\df Ul (x 0 , •) — d/u 2 (x o, -)IU p (o,t), (2.5) 

\\ dtU \( xQ , •) — d/u 2 (x o, -)||lp(o,t) < C||/i — /2|U°°(o,t)- (2.6) 

In particular, if we take p = oo in (2.2), then 

C 1 \\dfui(x 0 , •) — d/u 2 (x 0 , -)||l“(o,t) < ||/i — / 2 ||l°°(o,t) 

< C\\df Ul (x 0 , •) - dfu 2 (x 0 , -)||l-(o,t)- 

Theorem 2.2. Let condition (2.3) be fulfilled and v t be the solution of (1.2) for f — ./) G 


L°°(0,T) with ||/*|| l°°(o,t) < M (i — 1,2). We assume that there exist xq G D and 6 > 0 
such that 

\q(x 0 ,t)v 2 (x 0 ,t)\ >8, a.e. t G (0 ,T). (2.7) 

Then there exists a constant C > 0 depending on M, T, Q, 8 and ||g||L°°( 0i :r ; ,R 2 (O)) such that 

C 1 \\d?v 1 (xoi •) — dfv 2 (x 0 , OHz/^o.t) < ||/i — fi ||L°°(0,T) 
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< c\\d?vi(x 0 , •) - d/‘v 2 (x 0 , -)||l-(o,t). ( 2 . 8 ) 

In Theorem 4.4 of Sakamoto and Yamamoto [21], a similar problem to Theorem 2.1 is 
considered, but our result is more applicable in the point of view that the factor R(x, t ) is 
also allowed to depend on t. Moreover, we may assume less regularity for R in Theorem 
2.1. The arguments of Theorem 2.2 can also be applied to parabolic equations in order to 
consider the result of Theorem 1.1 in [7] with observations of the solution at a point xq G 
when C R d , d = 1, 2, 3. 

For such inverse problems with a = 1, we can also refer to Section 1.5 of Prilcpko, Orlovsky 
and Vasin [17], Cannon and Esteva [ 6 ] and Saitoh, Tuan and Yamamoto [19, 20], for example. 
In our main results, we assume conditions (2.4) and (2.7), which means that the observation 
point cannot be far from the source. On the other hand, in [6] and [19, 20], the determination 
of time dependent factor in the source term is considered without assuming such conditions 
and the logarithmic type and Holder type estimates are proved respectively. However, the 
results for fractional diffusion equations without these conditions have not been obtained 
yet. Here we restrict ourselves to the case with assumptions (2.4) and (2.7), and show the 
Lipschitz type stability. 

Let us remark that the results of this paper can be extended to the case d > 4. For this 
purpose additional conditions such as more regularity for a t j and dLl are required. In order 
to avoid technical difficulties, we only treat the case d < 3. 

3. Forward problem 

This section is devoted to the proof of unique existence and regularity of the solution of 
the IBVPs (1.1) and (1.2). 

Proposition 3.1. Let conditions (2.1) and (2.2) be fulfilled. Then the IBVP (1.1) admits 
a unique solution u G C([0, T]; H 2 (Ll)) satisfying 

Au G C([0, T\-H 2 \Q)) and d/u G L p ( 0, T; 
for all 0 < 7 < 1 — l/(pa). Moreover we have 

|^ m ||c([0,T , ];R' 2 RQ)) + \\dt U \\ LP^O.TiH 2 ^ (Q)) < C || fR || LP(0,T;H 2 (O)) • (3.1) 

with C > 0 depending on Ll, T and 7 


Proposition 3.2. Let conditions (2.1) and (2.3) be fulfilled. Then the IBVP (1.2) admits 
a unique solution v € C([0, T]; H 2 (Ll)) satisfying 

Av e C([0, T]; H 2 \LL)) and d?v e L°°( 0, T; 77 27 (H)) 
for all 0 < 7 < 1. Moreover, we have 

|“ / ^ ; ||c([0,T];iT 2 T'(f2)) + \\dt v \\ L °°(0, T ; H 2 'i(tt)) < C 11 Vo 11 H 4 (f2) 
withC depending onVt, T, ||/||z,°° ( 0 ,r), ||g||L->(o,T;^ 2 (^)) an &l- 


(3.2) 
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If all coefficients are independent of time variable t, then we can apply eigenfunction 
expansion and the problems can be reduced to ordinary differential equations of fractional 
order (e.g. [21]). However, since we consider the determination of the time dependent factor 
of coefficients, we apply fixed point theorem to show the unique existence of the solutions to 
(1.1) and (1.2) as in Beckers and Yamamoto [4]. 

In order to prove these results, we consider the IB VPs with more general data in the next 
subsections. 


3.1. Intermediate results. We introduce the following IB VPs 

d?u(x, t ) + Au(x, t) = F(x, t ), (x, t) G Q, 

B a u(x,t ) = 0, (x,t) G E, 

u(x, 0) = 0, iGfl, 

dfv(x, t ) + Av(x, t ) + p(x, t)v(x, t ) = F(x, t), (x, t ) G Q, 

B cr v(x,t ) = 0, (x,t) G E, 

v(x, 0) = 0, iGfl, 

and 

dfv(x, t) + Av(x, t) + p(x, t)v(x, t ) = 0, (x, t) G Q, 

B a v(x,t ) = 0, (x,t) G E, 

v(x, 0) = Vo(x)-, x G Q. 

We also consider the following conditions 

F G L p ( 0, T; H'\ fi)), - < p < oo and B a F = 0 on E 


a 


and 


(3,3) 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


1) p G L°°(0,T; H 2 (Q)) (and d u p — 0 on E if a ^ 0), 

2 ) vo G H 4 (Q) and B a v o = B a (Avo) = 0 on <90. 

Note that if we set F(x,t) = f(t)R(x,t), then conditions (2.1) and (2.2) are equivalent 
to (3.6). Similarly, if we assume p(x,t ) = f(t)q(x,t), then conditions (2.1) and (2.3) are 
equivalent to (3.7). Now let us consider the following intermediate results. 


Lemma 3.3. Let condition (3.6) be fulfilled. Then the IBVP (3.3) admits a unique solution 
u G C([0, T]; H 2 (Ll)) satisfying 

Au G C([0, T]; H 2 \Ll)) and dfu G L p ( 0 , T; H^(Ll)) 
for all 0 < 7 < 1 — l/(pa). Moreover we have 

\\Au\\c([0 y T];H 2 ~<(a)) + \\dt u \\ L P(p, T -, H ^(p.)) < C\\F\\ L p{p, T -,H 2 (Q)) 
with C > 0 depending on Ll, T and 7 . 


(3.8) 
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Lemma 3.4. Let F G L°°(0, T ; H 2 (Ll)) satisfy B a F = 0 and condition 1) of (3.7) be fulfilled. 
Then the IBVP (3.4) admits a unique solution v G C([0,T]-, H 2 (I})) satisfying 

Av G C([0, T]; tf 27 (fi)) and d?v G L°°( 0, T; tf 27 (fi)) 

for all 0 < 7 < 1. Moreover, we have 

||^||c([ 0 ,T];// 2 T(n)) + 11 V 11 L°° (0,T;H 2 1 (Q)) < C|| || L°°( 0 ,T;// 2 (Q)) (3.9) 

with C depending on T, ||p||z,°°(o,T ; ,R 2 (n)) and 7 . 


Lemma 3.5. Let condition (3.7) be fulfilled. Then the IBVP (3.5) admits a unique solution 
v G C([0, T]; H 2 (Ll)) satisfying 

Av G C([0, T]; H 2 \LL)) mid dfv G L°°(0, T; tf 27 (fl)) 
for all 0 < 7 < 1. Moreover, we have 

||*4t>||e([o,T];H 2 7(n)) + ||9f?;|| Loo(0jT . ff 2 7 (n)) < C||w||/nyn) (3.10) 

with C depending on Lt, T , ||p||L o °(o,T ; R' 2 (a)) an d 7- 


From these three lemmata we deduce easily Propositions 3.1 and 3.2. 

3.2. Preliminary. We define the operator A as A + l in L 2 (Vt) equipped with the boundary 
condition B a h = 0; 

ID(A ) := {h G H 2 (Ll)- B a h = 0 on 
[Ah :=Ah + h, h G D(A). 

Then A is a selfadjoint and strictly positive operator with an orthonormal basis of eigen¬ 
functions (4> n ) n > 1 of finite order associated to an increasing sequence of eigenvalues (A n ) n >i. 
We can define the fractional power A 7 , 7 > 0, of A by 


(3.11) 


D(A 7 ) := Ih G L 2 (M)- ^A 27 |(/q0 n )| 2 < 00 \ , 

l n= 1 

OO 

A 7 h:=J2 X n( h i<t>n)<f> n , heD(A 7 ). 


(3.12) 


71=1 


Then D(A 7 ) is a Hilbert space with the norm || • ||d(a^) defined by ||h||_D( 7 D) : = ||^4 7 /i||. 
Since D(A) is continuously embedded into H 2 (Vt) with norm equivalence (see Theorem 5.4 
in Chapter 2 of [13] for example), we see by interpolation that 

D(A 7 ) C H 27 (LI), 

C 1 ||h||H 2 T'(Q) < II^IId(at') < C , ||/?.||h 2 7 (q), h G D{A 7 ) 


for 0 < 7 < 1 . 
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In order to prepare for the arguments used in this paper, we consider the following Cauchy 
problem in L 2 (fl); 

dfu{t) + Au(t) = F(t ), t 

u(0) = 0. 

We define the operator valued function {S A (t)}t> o by 

OO 

S A (t)h = ^(Mn)£«,i(- \ n t a )(j> n , 

n =1 

with E aJ 3 , a > 0,/3 G M, the Mittag-Leffler function given by 


e(o ,T), 


(3.13) 


h g l 2 (n), t > o, 


E a ,p(z) — y] 


k =o 


T(ak + j3)' 


Recall that S A {t) € W 1 , 1 (0, T; B(L 2 (Q))) (e.g. [4] and [21]). Moreover, similarly to Theorem 
2.2 in [21], for F G L°°(0, T; L 2 (Q)), problem (3.13) admits a unique solution given by 

u{t)= [ A- 1 S , A (t-s)F(s)ds. (3.14) 

Jo 

This solution is lying in L°°(0, T; D^A 1 )) for 0 < 7 < 1, and, in view of Theorem 1 in [4], 
we have 

WA^S'^hW < C'f Q( 1 - 7 ) - 1 ||h||, h G L 2 (Q), t > 0. (3.15) 

In particular, the mapping t ha A~ l S' A {t) belongs to L q ( 0, T; £>(L 2 (fl))) if q G (1, cxd) satisfy 
q( 1 — a) < 1. Now we apply the following Young’s inequality to (3.14); 


Lemma 3.6. Let f G L p ( 0, T) and g G L q ( 0, T ) with 1 < p, q < 00 and l/p + l/q 
the function f * g defined by 

f*9(t ) := / f(t~ s)g(s)ds 
Jo 


belongs to C[0, T\ and satisfies 


1. Then 


\f*g(t)\ < ||/||L P (o,q||c/||L«(o,q, te [0, T\. 


Proof. Let / and g be defined by 


m ■= 



t G (0,T), 
t ^ ( 0 ,T), 


and g(t) 


g(t), t G (0, T), 
0, t£(0,T). 


Then applying the Young’s inequality for functions on M (see Exercise 4.30 in Brezis [5] or 
Appendix A in Stein [23] for example), we obtain the desired result. □ 
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Let p G (l,oo] be as in (3.6). Noting that A and A 1 S' A (t) commute, we see that for 
FeIA(0,T;D(A)), 

Au{t) = f A~ l S' A (t — s)AF(s)ds. 

Jo 

By p > 1/a and (3.15), the mapping t i-A- A^S'^t) belongs to L q (0, T; £>(L 2 (fl))) where 
q G [1, oo) satisfies 1/p + 1/q = 1. Therefore by Lemma 3.6, u belongs to C([0, T]; D(A)) 
and satisfies 

\\Au{t)\\< f ||dl" 1 ^(t-s)||||^lF(s)||ds < C [ (f-s) a - 1 ||F(s)|| D(A) ds (3.16) 


“ C (/ s(a ” 1) ^ S ) 7 l|F||iP (°’^» ^ Cta ~ 1 /P \\F\\LP { o,T- D (A))- 
Thus we can define the map FL : L p ( 0, T; D(A)) — > C([ 0, T]; D(A)) by 

FL(w)(t) := f A~ l S A (t — s)w(s)ds, w G 1/(0, T; D(A)). 


(3.17) 


(3.18) 


(3.19) 


By using these estimates, we will show the unique existence of the solution applying the 
fixed point theorem. 

3.3. Proof of Lemmata 3.3-3.5. 

Proof of Lemma 3.3. Let A be the operator defined by (3.11), then the IB VP (3.3) can 
be rewritten as 

<9 f Q w(f) + Au(t) = u(t) + F(t), t G (0, T), 

k u(0) = 0, 

where u(t) := u(-,t) and F(t) := F(-,t). Noting that F G 1/(0, T; D(A)) by (3.6), we see 
from (3.14) that the solution u of (3.19) satisfies 

u(t) = H(u)(t) + H(F)(t), t G (0,T), 

where the map FL is defined by (3.18). Therefore we will look for a fixed point of the map 
g : C([0, T\-D(A )) -A C([0, T]; D(A)) defined by 

G(w)(t) := H(w)(t) + H(F)(t ), w G C([0, T\\D{A)), t G (0, T). (3.20) 

By (3.16), for w G C([0, T]; D(A)), we have 

\\ / H(w)(t)\\ D{A) <C J (t - s) Q_1 ||w(s)|| D(A )ds < C (^J (f - s) Q_1 ds^ |MIc([o,t] ; d(A)) 


Ct c 


a 


Mlc([0,T];D(A))- 


Repeating the similar calculation, we get 

rt 


\\FL 2 w{t)\\ DlyA) <C J ( t-s) a 1 \\'Hw(s)\\ D ( A )ds < ^ {^J (t - s) a 1 s a ds^j |M|c([o,t] ; d(a)) 
C(T(a)t a ' 2 


T(2a + 1) 


| w ||c([0,T];D(A))- 
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By induction, we have 


\\ / H n w[t)\\ D {A) < 1 ll w llc([o,T];D(A))) w £ C([0, T\; D(A)). (3.21) 


T{na + 1) 


Therefore we obtain 


|| G H (wi) ~ G n (w2)\\c([0,T]-,D(A)) — \\'H n [w\ — W2 )\\c([0,T]-D(A)) 

C (T(a)T a ) n 

- “FT- Wl -W 2 C([0,T];D(A)) 

1 (not + 1) 

for wi,w 2 £ C([0,T]; D(A)). Since Q n is a contraction for sufficiently large n £ N, Q admits 
a unique fixed point u £ C([0, T]; D(A )) C C([0, T]; H 2 (Q)). Moreover we have 

n 

u = g( u ) = g n (u) = n n {u) + ^ n\F) 

k= 1 

for any n £ N. Now we estimate each 'H k (F). First, by (3.17), we have 

I|W(UWIId</i) < cr- l i-\\F\\ Lrlfs , T . MA)) . 

Next we apply (3.16) to have 

\\n 2 (F)(t)\\ DiA) <C fit-sy-'WUiFWWnwds 

Jo 

< C Qf (t - S )«-V-VP^ \\F\\lp( 0 ,T;D(A)) 

= r( 2 a+ 1 - 1 /p) * l|F|l " l "- TiD( ^) 


r( 2 a + 1 - 1 /p) 


Repeating the similar calculation, 


\mF)(t)\\DW < c / (t-sr-yn 2 (F)(s)\\ DiA) ds 


< 


CT(a) 


r( 2 a + 1 — 1 /p) Vi 0 
cr(q) 2 t 3a - 1 /p 
r(3a + 1 — 1/p) 


(■ t-s) a 1 s 2a 1/p ds ) IIF|| iP(0 , T;D(A)) 
||-F||lp(0,T;D(A))- 


By induction, we obtain 

||"H^(F) ||c([0,T];D(A)) < 

Therefore 


*, mll „ cr(a)‘- 1 r‘«- 1 /» 


T(ka + 1 — 1/p) 

n 

u\\c([0,T]-,D(A)) < ||'H n (u)||c([0,T];D(A)) + ^ \\F. k (F) ||c([0,T];D(A)) 


||-F||lp(0,T;D(A))- 


k=\ 
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„ C (r(a)T“)” „ , v' Crtaf-'T'™- 1 /? „ M1 

- r(»a +1) IMIc([w1iD(Jl » + a, r(to +1-1/p) l|F|ll ’ (0 ' TiDM)l 


and by taking sufficiently large n G N, we obtain 

IM|c([ 0 ,T];D(A)) < ^||-F1 |lp(0 ,T;.D(A)) 
with C depending on T and hi. 

Now fix 0 < 7 < 1 — 1 /(pa). Then for all t G (0, T), we have Au(t) G _D(H 7 ) with 


(3.22) 


A 1 (Au){t) = / A 7 ^^(t — s)(.Au(s) + AF(s))ds 
Jo 

and by (3.15), we have 

l^ 7 _ 1 ^(t)||B ( L a( n))<C'^ 1 , 

where // := a(l— 7 ). Since p > 1/p, the mapping t ha H 7-1 ^^) beongs to L 9 (0, T; £>(L 2 (fi)) 
where q G [1, 00 ) satisfies 1/p + 1/q = 1. Therefore Au belongs to C([0, T]; F>(H 7 )) and 

rt 


\\Au{t)\\ D ( A i) - 


A 1 - S A (t - s ) (aLu(s) + AF(s)) ds 


< C f (t-s) M 1 ||'u(s)|| D(A )ds + C' [ ( t-sY l 

Jo Jo 

< C y\t - sr^ds^j ||n|| C ( [0l T];D(A)) + C (j* sQ{fl ~ 1)ds ) 7 \\nLP( 0 ,t;D(A)) 

< CT^\\u\\c({o,t]-,d(A)) + CT P 1 //p ||T 1 ||lp(o,T;D(A)) (3.23) 

Combining this with (3.22), we have 

\\Au(t)\\ D (A-Y) < C\\F\\lp( 0,T]D(A)) < C , ||-P 1 ||z,P( 0 ,T;l/ 2 (n))- 
Hence we deduce that .Aw G C([0, T]; F 27 (fl)) and 

||*4w||e([0,T];.ff 2 7(n)) < ^11 II U p (0,T ;H 2 (Q)) • 

By the original equation dfu = — Au + F, we see that d/u belongs to L p ( 0, T ; dd 27 (fi)) with 
the estimate; 

||<9fu||z / p(0,T;H 2 7(Q)) < C||.Aw||i>(o,T;i/27(n)) + C\\F\\Lp( 0 ,T;H 2 ~i(£ 1 )) 

< C||.Aw||c([0,T];,H' 2 7(n)) + C , ||-f 1 ||LP( 0 ,T;fl’ 2 (n)) 

< C\\F\\lp(0,T-,H 2 (Q)), 

which implies (3.8). Thus we have completed the proof. □ 


For the proof of Lemma 3.4, we prepare the following fact; 


Lemma 3.7. Let u,v G H 2 (Q) and d < 3, then uv G H 2 (Q) with the estimate 

1 1 wu 1 1 # 2 (n) < C'||u||^2(q) 


with C depending on ||w||_ff 2 m). 



DETERMINATION OF COEFFICIENTS IN FRACTIONAL DIFFUSION EQUATIONS 


11 


For this lemma, see Theorem 2.1 in Chapter II of Strichartz [24], 


(3.24) 


Proof of Lemma 3.4. Similarly to Lemma 3.3, the IBVP (3.4) can be rewritten as 

d?v(t) + Av(t) = (1 - p(t))v(t) + F(t), 

^(0) = 0, 

where v(t) := v(-,t) and F(t) := Moreover p(t) denotes the multiplication operator 

by p(x,t). Then we can see that the solution v of (3.24) is a fixed point of the map /C : 
C([0, T]; D(A)) -> C([0, T]; D(A)) defined by 

K(w)(t):=(H{l-p{t))w)(t) + H(F){t), weC([0,T]-D(A)), te(0,T). 

Indeed, Lemma 3.7 and condition 1) of (3.7) yields that (1 —p)w belongs to L°°(0, T; D(A)) 
and satisfies 

11(1 -p(f)M*)||r>(A) < C\\w(t)\\ D ( A ) 

with C depending on \\p\\l°°(o,T;H 2 (q.))- Therefore we can see that /C maps C([0, T]; D(A)). 
Moreover, by the similar calculation to (3.21), we have 

C(T(a)t c 


11(77(1 -p))' l (^)l|c([0,T];D(A)) < 


Y{na + 1) 


|w||c([o,t];D(A))) w G C([0, T];D(A)) (3.25) 


and 


||(«(i -p))”- 1 («Ullc( [ o,T| i D W ) < r(n!r+“i ) )“ " F ll L °°< 0 . T;D (' < ) )' F € W’TMA))- 


By (3.25), we find 


(3.26) 


C(T(a)T a ) n 

\\IC n (wi) — /C n (w 2 )||c([ 0 ,T];D(A)) < ~pT- —11^1 — w 2 ||c([ 0 ,T];D(A)), 


T(na + 1 ) 

Wi,w 2 e C([0,T]; D(A))i 

which implies that 1C admits a unique fixed point v G C([0, T]; D(A)) C C([0, T]; H 2 (Q)). 
Then we have 

n 

v = K(v) = = (M(l -p(i)))» + £(W(1 -pm k -\nF)- (3.27) 

k =1 

Repeating the argument in the proof of Lemma 3.3, we deduce from (3.25), (3.26) and (3.27) 
that 

IMIc([o,t] ; d(A)) < C , ||-^ 1 ||l oc ’(o,T;D(A)) (3.28) 

with C depending on T, and \\p\\l°°(o,T;H 2 {Ci))- 
Next we fix 0 < 7 < 1. Similarly to (3.23), we have 

Av(t) e D(A^), t g (0,T) 

and 

rt 


\\Av(t)\\ D(A ^ < C 


A 7 S A (t — s) ((. 4.(1 — p(s))v)(s) + AF(s)) ds 
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<c [ ( t- s y i 1 (11(1 -p(s)Ms)||d(A) + II^(s)I|d ( A)) ds 
Jo 

<C f (t^s)^ 1 (|Ks)|| d(a) + ||i ;, (s)||£, ( A ) )d« 

Jo 

with n = a(l — 7 ). Therefore Av belongs to C([ 0, T]; 77 27 (f2)) and satisfies 

||v^'i ; ||c([0,T];_ff 2 7(a)) < ||v4n||c([0,T];D(A7)) < CT^ (||n||c([0,T];D(A)) + ||-F||l°°(0,T;D(A))) 

< C\\F\\l<x>(P,T-,D(A)) < ^II^I|l°°(0,T;JT 2 (£T)), 

where we have used (3.28). Moreover, combining this with the original equation, we also 
have dfv G L°°(0, T; 77 27 (f2)) and (3.9). □ 


Proof of Lemma 3.5. We split the solution v of (3.5) into two terms v = w +17 where w 
solves 

{ d^w(x, t ) + Aw(x, t ) + p(x, t)w(x, t ) = F(x, £), (x, f) G Q, 

B a w(x,t) = 0 , (x,t) e S, (3.29) 

w(x,0) = 0, 

with F(x,t ) := — (M + p(x, t))v 0 (x). Then (3.7) implies F G L°°(0, T; D(^4)) with the 
estimate 


||-^||l°°(0,T;// 2 (O)) < C\\v 0 \\h^(Q,)- 

By Lemma 3.4, the IBVP (3.29) admits a unique solution w G C([0, T]; 77 2 (f2)) satisfying 
Mw G C([0, T]; i/ 27 (f])) and <9> G L°°(0, T; i/ 27 (fl)). 


Moreover 


\\^ w \\c([0,T]-,H 2 'y(fl)) + \\dt W \\L°°(0,T-,H 2 ^(yi)) < C11 11 L°° (0,T-,H 2 (fl)) < C\\ v 0 || ■ 

Therefore the IBVP (3.5) admits a unique solution v G C([0, T]; H 2 (Q)) satisfying 
Av G C([0, T]; i/ 27 (f2)) and G L°°(0, T; iV 27 (0)). 

From the above estimate, we deduce (3.10). □ 


4. Proof of Theorem 2.1 

In this section, we prove Theorem 2.1. To this end, we prepare the following lemmata 
with Gronwall type inequalities; 


Lemma 4.1. Let C,a > 0 and u, d G L 1 (0, T ) be nonnegative functions satisfying 
u(t)<Cd(t) + C [ (t — s)“ _1 u(s)ds, i G (0,7), 


then we have 


u(t ) < Cd{t) + C (t- s) a ~ l d(s)ds, t G (0, T ) 
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For the proof, see Lemma 7.1.1 p.188 of [12]. 


Lemma 4.2. We take 2 < p < oo and p > 2/p. Let f G L°°(0,T) and u, R G L p (0,T) be 
non-negative functions satisfying the integral inequality 

f(t) < u(t) + [ (t — s) tl ~ 1 f(s)R(s)ds, a.e. t G (0,T). (4.1) 


Then we have 

||/||lp(0,T) < C'||m||lp(0 1 T), 

where the constant C depends on p, p, T and ||-R||lp(o,t)- 


(4.2) 


Proof. We set d(t) := ||/||x, P ( 0t )- From equation (4.1), we have 


which implies 


\f(s)\ p <C\u(s)\ p + C 


dft)<C\\u\\ p LP{0T) + C 


(s - 




ds. 


(4.3) 


Now we estimate the right-hand side of the above. By the Cauchy-Schwarz inequality, 

/ s r-s / rs \ l/ 2 / rs \ l/ 2 

|/«)UOr' 2 ^ = f l/K)l p/2 ' |«K)| P/ V{ < (j U |flK)|V{) , 

that is, 

\\fR\\LW H o,s) < ||/|Up(0,s)||-R||lp(0,s)- 

Therefore if p >2, then Lemma 3.6 yields that 

l/r 


(s - ar^nomdn 


<( / C^ds) ||/i?|| LP / 2(0 , s ) <C\\f\\ LP{0 , s) ||i?||LP(o, s) , 
0 / 


where r G [1, oo) satisfies 2/p + l/r = 1. For p = 2, we also have 

’(a -| L1(0 , S) < qiziu^ii^ii^^). 


Thus for any p > 2, we have 




< Cd(s), 


(4.4) 


where C depends on T, p, /i and ||-R||lp(o,t)- By (4.3) and (4.4), we have 

d(t) < C||-u||^ P (o T ) + C f d(s)ds, te(0,T). 

Jo 

Hence by the Gronwall inequality, we have 

d{t) < C\\u\\ p LP{QT)r t G (0, T) 

with C depending on p, p, T and ||-R||lp(o,t)- Thus we have proved (4.2). 


□ 



14 


KENICHI FUJISHIRO AND YAVAR KIAN 


Now we are ready to prove Theorem 2.1. 

Proof of Theorem 2.1. Let w* be the solutions to (1.1) corresponding to (i = 1,2) and 
set u, U\ — U ‘2 and / := f\ — f 2 . Then u solves (1.1) and is given by 

u(t) = f A~ l S' A (t — s)u(s) + f A~ l S' A (t — s)f(s)R(s)ds, 

Jo Jo 

where u(t) := u(-,t ) and R(t) := R(-,t). 

First we estimate ||w(t)||u(A)- Similarly to the calculation in (3.16), we have 

\\u(t)\\ D{A) < C j (t - s) Q_ 1 |Hs)||D(A)ds + C [ (t - s)“" 1 |/(s)|||i?(s)|| I )( J 4 )ds 

Jo Jo 

= c [ (t - s) Q_1 ||M(s)||DM)ds + edit), (4.5) 


where we have set 

d{t) := [ (t-s)“^ 1 |/(s)||| J R(s)||D(A)ds. 

Jo 

Applying Lemma 4.1 to (4.5), we have 

\\u(t)\\D(A) < Cd{t) + C f (t — s) Q_ 1 d(s)ds, 0 < t < T. 


Here for v > 0, we note 


(t — s) u 1 d(s)ds = / (t — s) L 


(s- o a -vmm)\\D { A)dt; ds 


(t-sy-^s-^ds |/(0llltf(0IIW£ 


'o ve 


= B(u,a) / (4.7) 

Jo 

where B(-, •) is the Beta function. In particular, for v = a, we have 

f it-s) a ~ 1 d{s)ds = B(a,a) [ (t - s) 2a_1 |/(s)|||R(s)||r>(A)ds 
Jo Jo 

<T a B(a,a) Ai - s) a - 1 |/(s)|||i2(s)|| D(A )ds 
Jo 

< Cd{t). 

Hence the following estimate follows from (4.6); 

||w(f)||r»(A) < Cd(t), 0 <t<T. (4.8) 

Next we estimate ||.4 /u(£)||.d(ay) for d/4 < 7 < 1 — 2 /(pa). Repeating the calculation in 
(3.23), we find 

\\Au{t)\\ D{A ^ < C [ (t- s)^ 1 (||«(s)||£, ( a) + |/(s)|||#(s)||d ( a)) ds, a.e. t G (0, T) 
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where /i = a( 1 — 7 ). By (4.7) with v = /i and (4.8), we obtain 

\\Au(t)\\ D ( A y) < C f (t - sY~ 1 d(s)ds + C [ (t - s) M_ 1 |/(s)|||i?(s)|| D ( A )ds 
Jo Jo 

= CB(fi,a) f\t - s) |,+ “- 1 |/(s)|||if(s)ll 
Jo 

+ C fit-sY-^f^lWR^Wn^ds 
Jo 

< CT a B(n, a) [ (t-sr-'lfisMRisYn^ds 
Jo 

+ C [ (t-sr-'lfisMRYYn^ds 
Jo 

<C ['(t-8)^\f(a)\\\R(a)\\ DiA) d8. 

Jo 

Finally we estimate \Au{xo,t)\ and complete the proof. Since 7 > d/ 4, the Sobolev 
embedding theorem yields 

\Au(x 0 ,t)\ < C\\Au(-,t)\\ H 2 y {n) < C\\Au{t)\\ D{Al) < C [ (t- s)^ 1 |/(s)|||i?(s)|| D(A) ds. 

Jo 

(4,9) 


From the original equation, we get 

f(t)R(x 0 ,t ) = d^u(x 0 ,t) + Au(x 0 ,t), a.e. t G (0,T). (4-10) 

Combining this with (2.4) and (4.9), we get 

\f(t)\ < ^ (\^u(x 0 ,t)\ + \ Au(x 0 ,t)\) 

<C\d?u(xo,t)\ + C fit-sy-^fisMRtfWDwds, a.e. t G (0, T) (4.11) 
Jo 

with C depending on 6 , 11 and T. By Lemma 4.2, we see that 

||/||lp(o,t) < C\\d?u(x 0 , -)||lp(o,t), 
which implies (2.5). Moreover, by (4.9) and (4.10), we have 
\d?u(x 0 ,t)\ < \f(t)R(x 0l t)\ + \ Au(x 0 ,t)\ 

<C\fm\R(;t)\\m(n) + C [\t - 3 )^( 3 )\\\R(s)\\ D{A) ds 

Jo 

< C , ||/||l°°(0,T)||-R(^)||d(A) + C'||/|| L oo(0,D [ (t — sY 1 ||l?(-s)||D(A)ds 

Jo 

< C\\f\\L°°(o,T)\\R(t)\\D(A) + (7||/||l°°(o,t)^ 1 ^ p \\R\\lp(o,T;D(A)), a.e. t G (0,T). 


Therefore, 

\\dfU(x 0 , -)||lp( 0 ,T) < (7||/||l°°(0,T)||7?||lp(o,T;D(A)) + C'||/|| L oo( 0iT )T M ||i?|| L p(o i T;D(A)) 
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< C\\f\\L°°(0,T)- 

Thus we have proved (2.6). □ 

5. Proof of Theorem 2.2 

In this section, we prove Theorem 2.2. We first prepare the following generalized Gron- 
wall’s inequality; 

Lemma 5.1. Let p, a, b > 0 and f G L 1 (0, T ) be nonnegative function satisfying the integral 
inequality 

f(t)<a + b f ft — s) M_1 f(s)ds, a.e. t G (0, T). 

Jo 

Then we have 

fit) < aE M ((ftTM) 1 /^) , a.e. t G (0,T). 

For the proof, see Lemma 7.1.2 on p.189 of [12]. Now we are ready to prove Theorem 2.2. 

Proof of Theorem 2.2. Let v l be the solutions to (1.2) corresponding to fi {i = 1,2) and 
set v := V\ — v -2 and / := f 2 — fi- Then v solves (3.4) with p(x,t ) = fi(t)q(x,t) and 
F(x,t ) = f(t)q(x,t)v 2 (x,t). Recall that v is given by 

v(t) = f A -1 S' A {t - s){{ 1 - p{t))v)(s) + I f(s)A- 1 S' A (t-s)R(s)ds, 

Jo Jo 

where we have set v(t) := v(-,t) and R{t) := q(-, t)v 2 (-, t). Moreover, p(t) denotes the 
multiplication operator by p(x,t) := fi(t)q(x,t). 

First we estimate ||f(^)IId(A)- Since (1 — p(t))v(t), R(t) G D(A) by (2.3), we repeat the 
calculation in (3.23) to have 

||F(t)|| D(A) < C [ (t-s) Q_ 1 ||(l -p{t))v(s)\\ D (A)ds + C [ it - s)"" 1 |/(s)|||i?(s)|| i? ( J 4)ds 
Jo Jo 

< c [ it - s) q- 1 ||u(s)|| d( a)gG + C [ it - s) Q- 1 |/(s)|ds. 

Jo Jo 

with C depending on Q, M and ||q , ||L°o(o,T;H 2 (n))- Then repeating the arguments used in 
Theorem 2.1, we obtain 

(t — s)“ _ 1 |/(s)|ds, 0 < t < T. 

and from this estimate we also deduce that for any 0 < 7 < 1 , 

\\Av{t)\\ D{A y) < C [ it - s) M_ 1 |/(s)|ds, 0 < f < T, 

Jo 

where p, o(l — 7 ). Therefore by taking 7 G (d/4,1), we have 

I Av(x 0 ,t) +pix 0 ,t)vix 0 ,t )| < C\\Av(-,t) +p(-,t)v{-,t)\\ H ^(n) 

<C\\Avi;t )\\ H2 

na) + C\\vi-,t )\\ H 2 7(0) 
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<C\\Av(t)\\ D ^<C (t-sY 1 |/(s)|ds. (5.1) 

Jo 

From the original equation, we have 

f(t)R(x 0 , t ) = d?v(x 0 , t ) + Av(x o, t ) + p(x 0 , f)u(rc 0 , £), a.e. f G (0, T). (5.2) 

On the other hand, from (2.7), we deduce that 

|R(x 0 , £)| > c > 0, a.e. t G (0, T) 

with c depending on <5, O and T. Therefore, combining this with (5.1) and (5.2), we obtain 
\f(t)\ < C\d?v(x 0 ,t)\ +C\Av(x 0 ,t)+p(x 0 ,t)v(x 0 ,t)\ 

< C\\d?v(x 0 ,-)\\ L <~ { 0 ,T) + C [ (f-s) M_1 |/( s )Ms, a.e. t G (0, T). 

Jo 

Applying Lemma 5.1, we see that 

\f(t)\ < C\\d^v(x 0 , •)||l oo (o,t)- 

Thus we have proved the second inequality in (2.8). Moreover, by (5.2), we have 

\<%v{x Q ,t)| < \f(t)R(x 0 ,t)\ + | Av(x 0 ,t) +P(x 0 ,t)v(x 0 l t)\ 

<\f(t)\\\R(-,t)\\ Di A) + C [ (t- sY-'lfis^ds 

Jo 

T>'\ 


□ 


< C ( \\R\\l°°(o,t-,d(A)) H-] ||/||l°°(o,t)- 

I 1 J 


Thus we have proved the first inequality in (2.8). 
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